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Abstract. The Quasi-Elastic (QE) contribution of the nuclear inclusive electron model developed in reference 1 1] is extended 
to the study of neutrino/antineutrino Charged Current (CC) and Neutral Current (NC) induced nuclear reactions at intermediate 
energies. Long range nuclear (RPA) correlations, Final State Interaction (FSI) and Coulomb corrections are included within 
the model. RPA correlations are shown to play a crucial role in the whole range (100-500 MeV) of studied neutrino energies. 
Results for inclusive muon capture for different nuclei through the Periodic Table are also discussed. In addition, and by 
means of a Monte Carlo cascade method to account for the rescattering of the outgoing nucleon, we also study the CC and 
NC inclusive one nucleon knockout reactions off nuclei. 
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INTRODUCTION 

Neutrino physics is at the forefront of current theoretical 
and experimental research in astro, nuclear, and particle 
physics. Indeed, neutrino interactions offer unique op- 
portunities for exploring fundamental questions in these 
domains of the physics. One of these questions is the 
neutrino-oscillation phenomenon, for which there have 
been conclusive positive signals in the last years j^. 
The presence of neutrinos, being chargeless particles, can 
only be inferred by detecting the secondary particles they 
create when colliding and interacting with matter. Nuclei 
are often used as neutrino detectors, thus a trustable inter- 
pretation of neutrino data heavily relies on detailed and 
quantitative knowledge of the features of the neutrino- 
nucleus interaction. There is a general consensus among 
the theorists that a simple Fermi Gas (FG) model, widely 
used in the analysis of neutrino oscillation experiments, 
fails to provide a satisfactory description of the measured 
cross sections, and inclusion of further nuclear effects is 
needed Q3I] . 

Any model aiming at describing the interaction of neu- 
trinos with nuclei should be firstly tested against the ex- 
isting data on the interaction of real and virtual pho- 
tons with nuclei. For nuclear excitation energies rang- 
ing from about 100 MeV to 500 or 600 MeV, three dif- 
ferent contributions should be taken into account: i) QE 
processes, ii) pion production and two body processes 
from the QE region to that beyond the A( 1232) resonance 
peak, and iii) double pion production and higher nu- 
cleon resonance degrees of freedom induced processes. 
The model developed in Refs. [Q]] (inclusive electro- 
nuclear reactions) and [4] (inclusive photo-nuclear reac- 
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FIGURE 1. Diagrammatic representation of the neutrino 
selfenergy in nuclear matter. 



tions) has been successfully compared with data at inter- 
mediate energies and it systematically includes the three 
type of contributions mentioned above. Nuclear effects 
are computed starting from a local FG picture of the nu- 
cleus, an accurate approximation to deal with inclusive 
processes which explore the whole nuclear volume [4], 
and their main features, expansion parameter and all sort 
of constants are completely fixed from previous hadron- 
nucleus studies (pionic atoms, elastic and inelastic pion- 
nucleus reactions, A— hypernuclei, etc..) 05, 6]. Thus, 
and besides the photon coupling constants determined in 
the vacuum, the model of Refs. [ 1] and [4] has no free 
parameters, and thus the results presented in these two 
references are predictions deduced from the nuclear mi- 
croscopic framework developed in Refs. [5] and Jg]. In 
this talk, we extend the nuclear inclusive QE electron 
scattering model of Ref. 1 1 1. including the axial CC 10] 
and NC [8] degrees of freedom, to describe neutrino and 
antineutrino induced nuclear reactions in the QE region. 




FIGURE 2. Set of irreducible diagrams responsible for the polar- 
ization (RPA) effects in the Iplh contribution to the W— selfenergy. 
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FIGURE 3. V e — (top) and v e — (bottom) inclusive QE differential 
cross sections in oxygen as a function of the transferred energy, at two 
values of the transferred momentum. We show results for relativistic 
(REL) and non-relativistic nucleon kinematics. In this latter case, we 
present results with (FSI) and without (NOREL) FSI effects. For the 
three cases, we also show the effect of taking into account RPA corre- 
lations and Coulomb corrections (lower lines at the peak). 



We also present results for the QE (V/,V/iV), (Vi,l~N), 
(vi,ViN) and (V[,l + N) reactions in nuclei. We use a 
Monte Carlo (MC) simulation method to account for the 
rescattering of the outgoing nucleon. The first step is the 
gauge boson (W ± and Z° ) absorption in the nucleus. 
We take this reaction probability from the microscopical 
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FIGURE 4. Muon neutrino inclusive QE differential cross sec- 
tions in oxygen as a function of the transferred energy. The notation for 
the theoretical predictions is the same as in Fig. [J] 

many body framework developed in Refs. for CC 

and NC induced reactions. Some calculations found in 
the literature use the plane wave and distorted wave im- 
pulse approximations (PWIA and DWIA, respectively), 
including or not relativistic effects. The PWIA consti- 
tutes a poor approximation, since it neglects all types of 
interactions between the ejected nucleon and the residual 
nuclear system. The DWIA describes the ejected nucleon 
as a solution of the Dirac or Schrodinger equation with 
an optical potential obtained by fitting elastic proton- 
nucleus scattering data. The imaginary part accounts for 
the absorption into unobserved channels. This scheme is 
incorrect to study nucleon emission processes where the 
state of the final nucleus is totally unobserved, and thus 
all final nuclear configurations, either in the discrete or 
on the continuum, contribute. The distortion of the nu- 
cleon wave function by a complex optical potential re- 
moves all events where the nucleons collide inelastically 
with other nucleons. Thus, in DWIA calculations, the 
nucleons that interact inelastically are lost when in the 
physical process they simply come off the nucleus with 
a different energy, angle, and maybe charge, and they 
should definitely be taken into account. A clear example 
which illustrates the deficiencies of the DWIA models 
is the neutron emission process: (v/,/ _ n). Within the im- 
pulse approximation neutrinos only interact via CC inter- 
actions with neutrons and would emit protons, and there- 
fore the DWIA will predict zero cross sections for CC 
one neutron knock-out reactions. However, the primary 
protons interact strongly with the medium and collide 
with other nucleons which are also ejected. As a con- 
sequence there is a reduction of the flux of high energy 
protons but a large number of secondary nucleons, many 
of them neutrons, of lower energies appear. 

Finally, we would like to mention that the we have 
already started working on one and two pion production 
processes in nuclei, aiming to extend the present study 
beyond the QE peak to the A(1232) resonance region. 
As a first step, we have studied the production off the 



nucleon and some results have been also reported 

at this conference 



INCLUSIVE QE CROSS SECTIONS 

We will expose here the general formalism focusing on 
the neutrino CC reaction. The generalization to antineu- 
trino CC, neutrino and antineutrino NC reactions or in- 
clusive muon capture is straightforward. In the labora- 
tory frame, the differential cross section for the process 
V[(k) + A z ^r {k 1 ) +X reads: 
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with L and W the leptonic and hadronic tensors, respec- 
tively. On the other hand, the inclusive CC nuclear cross 
section is related to the imaginary part of the neutrino 
self-energy (see Fig. [TJ in the medium by: 



= -4- / lmL v (k:p(r))d 3 r 



(2) 



We obtain the imaginary part of the neutrino self-energy 
in the medium, Im£ v , by means of the Cutkosky's rules: 
in this case we cut with a vertical straight line (see Fig. [O 
the intermediate lepton state and those implied by the W- 
boson medium self-energy. Those states are placed on 
shell by taking the imaginary part of the propagator, self- 
energy, etc. We obtain for k° > 
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and thus, the hadronic tensor is basically an integral over 
the nuclear volume of the W— selfenergy (Yly^(q;p)) 
inside the nuclear medium. We can then take into ac- 
count the different in-medium effects and reaction mech- 
anism modes (W absorption by one nucleon or by a 
pair of nucleons, pion production, excitation of reso- 
nances,...) by including the correspondent diagrams in 
the W— selfenergy (shaded loop of Fig. [TJ. Further de- 
tails can be found in Refs. (HQ]. 

The virtual W gauge boson can be absorbed by one 
nucleon, lplh nuclear excitation, leading to the QE con- 
tribution to the nuclear response function. We consider a 
structure of the V — A type for the W + pn vertex, and use 
PCAC and invariance under G-parity to relate the pseu- 
doscalar form factor to the axial one and to discard a term 
of the form (/?^ + p'^)Y5 m the axial sector, respectively. 
Besides, and thanks to isospin symmetry, the vector form 
factors are related to the electromagnetic ones. We find 
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with the local Fermi momentum kf(r) = 
(37T 2 p(r)/2) 1//3 , M the nucleon mass, and E(p) = 
\/M 2 + p 2 . We work on an non-symmetric nuclear 
matter with different Fermi sea levels for protons, k F , 
than for neutrons, k n F (equation above, but replacing 
p /2 by p p or p„, with p = p p + p„). Finally, A^ v is 
the CC nucleon tensor (0] and is determined by the 
W + pn form-factors. The d?p integrations above can be 
done analytically and all of them are determined by the 
imaginary part of isospin asymmetric Lindhard function, 
U(q,k' F ,k F ). Explicit expressions can be found in |7[], 
where also antineutrino induced CC cross sections 
and the inclusive muon capture process in nuclei are 
discussed. Expressions for the NC hadron tensor can be 
found in Ref. ||8[]. 



Nuclear Corrections 

We take into account polarization effects by sub- 
stituting the particle-hole (lplh) response by an 
RPA response consisting of a series of ph and 
A-h excitations, as shown in Fig. [2] We use an 
effective Landau-Migdal ph-ph interaction 11211 : 
V = co{/o+/oTiT2+go5ia2+goCTi02TiT 2 }. In the 
vector-isovector channel (ggxt operator) we use an 
interaction |H 0, @1 with explicit n— meson (longitudi- 
nal) and p —meson (transverse) exchanges, and that also 
includes A(1232) degrees of freedom. This effective in- 
teraction is non-relativistic, and then for consistency we 
will neglect terms of order ff(p 2 /M 2 ) when summing 
up the RPA series. 

We also ensure the correct energy balance of the differ- 
ent studied processes by modifying the energy conserv- 
ing 8 function in Eq. (0]) to account for the experimental 
Q— value of the reaction. Besides, we consider the effect 
of the Coulomb field of the nucleus acting on the ejected 
charged lepton, as well. This is done by including the 
lepton self- energy Ec = 2k Vc(r) in the intermediate 
lepton propagator of Fig.Q] 

Finally, we take into account the modification of nu- 
cleon dispersion relation in the medium (FSI) by us- 
ing nucleon propagators properly dressed with a realistic 
self-energy. We use a dynamical model in which the nu- 
cleon self-energy depends explicitly on the energy and 
the momentum [13]. Thus, we compute the imaginary 
part of the Lindhard function (ph propagator) using real- 
istic spectral functions Spj,(co,p;p). 
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FIGURE 5. Ratio of inclusive QE cross sections induce by neu- 
trinos (top) and antineutrinos (bottom) in oxygen, as a function of the 
incoming neutrino energy. Solid and dashed lines stand for the predic- 
tions of the model of Ref. |7] and from the FG model, respectively. 



CC AND NC INCLUSIVE ONE 
NUCLEON KNOCKOUT REACTIONS 

We use a cascade method to account for the rescattering 
of the outgoing nucleon. The first step is the gauge bo- 
son (W* and Z° ) absorption in the nucleus, we take this 
reaction probability from the microscopical many body 
framework outlined in the previous section for CC and 
NC induced reactions. It is given by the inclusive QE 
cross sections d 2 o/dQ!dE' {Ql, E' are the solid angle 
and energy of the outgoing lepton) for a fixed incoming 
neutrino or antineutrino laboratory energy. We also com- 
pute differential cross sections with respect to d 3 r. Thus, 
we also know the point of the nucleus where the gauge 
boson was absorbed, and we can start from there our MC 
propagation of the ejected nucleon. After the absorption 
of the gauge boson, we follow the path of the ejected nu- 
cleon through its way out of the nucleus using a MC sim- 
ulation to account for the secondary collisions. Details on 
the MC simulation can be found in |8]. This MC simula- 
tor has been tested in different physical situations. It was 
first designed for single and multiple nucleon and pion 
emission reactions induced by pions ffl [TBll and has 
been successfully employed to describe inclusive (7, %), 



{y,N), (y,AW),..., {y,Nn) fiHH, {e,e'n), {e,e'N), 
(e,e'NN),..., (e,e'Nn),... d reactions in nuclei or the 
neutron and proton spectra from the decay of A hypernu- 
cleild. 

RESULTS 
Inclusive cross sections 

At low energies, the model provides a reasonable de- 
scription of the LSND measurement of the reaction 12 C 
(y t i,jX^)X near threshold, and of the accurate nuclear 
inclusive muon capture rates through the whole Peri- 
odic Table. Pauli blocking and the use of the correct 
energy balance provide the bulk of the corrections, but 
an accurate description of data is only achieved once 
the RPA correlations, including Ah excitations, are taken 
into account. Our approach provides one of the best ex- 
isting combined descriptions of the inclusive muon cap- 
ture in 12 C and the LSND measurement of the reaction 
12 C {Vjji,ll^)X near threshold [7], and certainly compa- 
rable to that achieved by other models, which implement 
a more sophisticated treatment of the dynamics of finite 
nuclei (see for instance the discussion in [20]). 

At intermediate energies the predictions of our model 
become reliable not only for integrated, but also for dif- 
ferential cross sections. We present results for incoming 
neutrino energies within the interval 150-500 MeV for 
electron and muon species (some results at higher ener- 
gies can be found in B2l|p . In Figs. [3] and |4] RPA and 
FSI effects on differential cross section are shown. RPA 
effects are extremely important in this range of energies 
and induce important corrections to the simple FG de- 
scription. On the other hand, FSI provides a broadening 
and a significant reduction of the strength of the QE peak. 
Nevertheless the integrated cross section is only slightly 
modified. Though FSI change importantly the shape of 
the differential cross sections, it plays a minor role for 
totally integrated cross sections. When medium polar- 
ization effects are not considered, FSI provides signif- 
icant reductions (15-30%) of the cross sections. How- 
ever, when RPA corrections are included the reductions 
becomes more moderate, always smaller than 7%, and 
even there exist some cases where FSI enhances the cross 
sections. This can be easily understood by looking at 
Fig. [4] There, we see that FSI increases the cross sec- 
tion at high energy transfers. But for nuclear excitation 
energies higher than those around the QE peak, the RPA 
corrections are certainly less important than in the peak 
region. Thus, the RPA suppression of the FSI distribution 
is significantly smaller than the RPA reduction of the dis- 
tribution determined by the ordinary Lindhard function. 

We have estimated the theoretical uncertainties of our 



model by MC propagating the uncertainties of its differ- 
ent inputs into differential and total cross sections ll22ll . 
We conclude that our approach provides QE v(v)- 
nucleus cross sections with relative errors of about 10- 
15%, while uncertainties affecting the ratios cj(jU) / c(e) 
and a{pL)/o(e) would be certainly smaller, not larger 
than about 5%, and mostly coming from deficiencies of 
the local FG picture of the nucleus [22]. Though nuclear 
corrections cancel importantly in these ratios, there still 
exist some effects ( Fig. [5]). 
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Nucleon Emission Reactions 

CC and NC nucleon emission processes play an im- 
portant role in the analysis of oscillation experiments. In 
particular, they constitute the unique signal for NC neu- 
trino driven reactions. Different distributions for both NC 
and CC processes can be found in [8], as example, we 
show here results for NC nucleon emission from argon 
(Fig. [6). The rescattering of the outgoing nucleon pro- 
duces a depletion of the high energy side of the spec- 
trum, but the scattered nucleons clearly enhance the low 
energy region. Our results compare well with those of 
Refs. I23I UM obtained by means of a transport model. 
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FIGURE 6. Neutral current 40 Ar(v,v+N) at 500 MeV 
cross sections as a function of the kinetic energy of the final 
nucleon. The dashed histograms show results without rescat- 
tering (PWIA) and the solid ones have been obtained from the 
MC cascade simulation. 
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